Quantum-jumps and photon-statistic in fluorescent systems coupled to classically 

fluctuating reservoirs 
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In this paper, we develop a quantum-jump approach for describing the photon-emission process of 
single fluorophore systems coupled to complex classically fluctuating reservoirs. The formalism relies 
on an open quantum system approach where the dynamic of the system and the reservoir fluctuations 
are described through a density matrix whose evolution is defined by a Lindblad rate equation. For 
each realization of the photon measurement processes it is possible to define a conditional system 
state (stochastic density matrix) whose evolution depends on both the photon detection events 
and the fluctuations between the configurational states of the reservoir. In contrast to standard 
fluorescent systems the photon-to-photon emission process is not a renewal one, being defined by 
a (stochastic) waiting time distribution that in each recording event parametrically depends on 
the conditional state. The formalism allows calculating experimental observables such as the full 
hierarchy of joint probabilities associated to the time intervals between consecutive photon recording 
events. These results provide a powerful basis for characterizing different situations arising in 
single-molecule spectroscopy, such as spectral fluctuations, lifetime fluctuations, and light assisted 
processes. 

PACS numbers: 42.50. Lc, 42.50.Ct, 42.50.Ar, 33.80.-b 



I. INTRODUCTION 

A powerful theoretical formalism called the quantum- 
jump approach |ll4ll| was introduced by the quantum op- 
tics community for describing experimental realizations 
of single open quantum systems subjected to a continu- 
ous measurement process. Even when only one system 
is under observation, the quantum-jump approach allows 
to define a system state (wave vector or density matrix) , 
whose dynamic takes into account our change of informa- 
tion due to the continuous measurement action. Apart 
from new insights in the quantum measurement theory, 
the quantum-jump approach provides an alternative for- 
malism for characterizing the radiation pattern of single 
fluorescent systems driven by a laser field. 

While a wide class of quantum optical systems can be 
studied with the quantum-jump approach [l|-y|i it has 
been scarcely applied in the context of single-molecule 
(fluorescence) spectroscopy (SMS) p^ - [l4| . i.e., in the 
characterization of single fluorescent systems coupled to 
complex host classically fluctuating environments, such 
as of those associated to biological or artificially designed 
nanoscopic reservoirs. The main task of SMS is to deduce 
the underlying environment stochastic dynamic from the 
statistical properties of the scattered laser field p^| - |25| . 
In most of the experiments, the scattered electromagnetic 
field is measured with photon detectors. Hence, it can be 
resolved photon-to-photon. 

For direct photon-detection measurement schemes the 
quantum-jump approach associate to each photon record- 
ing event a sudden disruptive change (wave vector col- 
lapse) in the system state, while in the middle intervals 
between consecutive events the-conditional-system evo- 
lution is smooth and non-unitary [l|-l6|. The formalism 



provides a simple technique for calculating and reproduc- 
ing the photon recording process. For Markovian dissi- 
pative dynamics the emission process is a renewal one, 
i.e., the statistic of the (random) time intervals between 
consecutive photon emissions is always the same, being 
defined by a probability distribution called waiting time 
distribution [l|, |5| . 

The main obstacle for applying the quantum-jump ap- 
proach for modeling SMS experiments comes from the de- 
scription of the environment fluctuations. As in general 
a full microscopic description is lacking, the complexity 
of the environment is taken into account by introducing 
effective time-dependent stochastic variables that may 
modify (parametrize) both the unitary and dissipative 
fluorophore evolution. In the context of the quantum- 
jump approach, it is not clear how these extra (classical) 
fluctuations must to be introduced or consistently inter- 
preted in terms of a continuous measurement action. 

On the basis of stochastic models, the formalism of 
generalized Bloch equations [26l - (30 | allows to determin- 
ing the photon counting probabilities, i.e., the probabil- 
ities of detecting n-photons up to a given time. Never- 
theless, from that approach it is not easy to know how 
the renewal property is broken by the external fluctua- 
tions, neither is known which kind of stochastic dynamic 
may reproduce the photon emission process. Then, ob- 
jects like the hierarchy of joint probabilities associated to 
the time intervals between consecutive photon detections 
events is also unknown. These statistical objects can be 
obtained, for example, from a time average along a single 
measurement trajectory [see Eqs. (fTU]) and (|lip ]. 

The main goal of this paper is to demonstrate that 
SMS experiments can be consistently described in the 
context of a quantum-jump approach. A general for- 



malism that allows to characterize the photon-to-photon 
emission process for a broad class of environment fluctu- 
ations arising in SMS is developed. In each case, we pro- 
vide (non-renewal) stochastic processes that reproduce 
the statistic of the photon recording events. The average 
of their associated dynamic in the system-bath Hilbert 
space recover the density matrix evolution. As a central 
result, we get explicit analytical expressions for the set 
of joint probabilities densities defining the statistics of 
the time intervals between successive photon recording 
events. Therefore, our analysis allows to quantify how 
and how much the photon emission process departs from 
a renewal one. 

The formulation of an alternative description of SMS 
experiments based on a quantum-jump approach relies 
on the possibility of describing both the fluorophore and 
the environment fluctuations through a density matrix 
formalism. In Ref. [3l| it was demonstrated that a broad 
class of SMS experiments can be studied through an open 
quantum system approach. The density matrix evolution 
is given by a Lindblad rate equation |32| , which allows to 
characterize in a unified way both the quantum nature 
of the fluorescent system as well as the classical nature 
of the environment fluctuations. Based on those results, 
which are consistent Isil with the formalism of stochastic 
Bloch equations [26l430l |. we formulate the present treat- 
ment. 

We remark that a similar analysis was developed in 
Ref. [33|. In contrast, our present analysis allows get- 
ting explicit expressions for the photon emission statis- 
tic, which is also analyzed in the limit of slow and fast 
environment fluctuations. Furthermore, an explicit for- 
mulation of the underlying stochastic photon emission 
process is presented. On the other hand, our results also 
clarifies some of the assumptions introduced in previous 
author's works |34]4j19 as well as in other stretched re- 
lated contributions [37]. 

The paper is outlined as follows. In Sec. II, on the basis 
of the results developed in Ref. [3l[ , we define the under- 
lying density matrix formulation. In Sec. Ill, we develop 
the quantum-jump approach. Both, the stochastic dy- 
namic and the statistical characterization of the photon 
emission process are established. In Sec. IV we apply the 
formalism for the case in which the measurement appa- 
ratus only gives information about the photon emission 
events. Different specific cases, such as lifetime fluctua- 
tions and light assisted processes, are analyzed in detail. 
In Appendix A we analyze the case of measurements that 
provide information of both the photon recording events 
and about the configurational reservoir transitions. In 
Sec. V we provide the conclusions. 



II. DENSITY MATRIX EVOLUTION 

The description of SMS experiments based on a den- 
sity matrix formalism relies on the possibility of finding 
analytically manageable microscopic interactions able to 



describe the environment fluctuations as well as their dy- 
namical influence over the system. In Ref. 1 3 111 , following 
an argument developed by van Kampen [38|, we mod- 
eled the environment through a set of (effective, coarse 
grained) macrostates, each one representing the manifold 
of quantum bath states that lead to the same system dy- 
namic. Then, the total microscopic dynamic is written in 
an effective Hilbert space defined by the external prod- 
uct of the Hilbert spaces of the system, the background 
electromagnetic field, and the configurational space asso- 
ciated to the bath macrostates. The system is modeled 
by a two-level optical transition whose characteristic pa- 
rameters, i.e., transition frequency and electric dipole, 
depend on the state of the environment. The dielectric 
constant of its local environment also is parametrized by 
the bath macrostates. After tracing out the electromag- 
netic field and the configurational states, the density ma- 
trix ps{t) of the system can be written as |31j 



psit) = T.Z7p^^^^ 



(1) 



Each auxiliary state pR{t) define the system dynamic 
given that the reservoir is in the _R-configurational bath 
state. i?max is the number of configurational states. The 
probability -Pfl(t) that the environment is in a given state 
at time t follows from 



PR[t)^TYs[pR{t)l 



(2) 



where Tr5[- • •] denotes a trace operation in the system 
Hilbert space. Therefore, the set of states {pR{t)} encode 
both the system dynamic and the fluctuations of the en- 
vironment. Their dynamic is defined by a Lindblad rate 
equation [S^l 

^^ = ^[HR,PR{t)]-^R{{D,PR{t))+-J[pRm{2,) 
- $]^{AtA,p;,,(t)}++ Y.^RR,,ApR,{t)A\ 



R' 
R'^R 



R' 
R'^R 



The first line of this equation define the unitary and dis- 
sipative system dynamic given that the bath is in the 
configurational state R. The Hamiltonian Hr reads 



Hr^ 



hujR, 



hn, 



2 a, + ^(ate— -* + ae+--*), (4) 



where 



^R = (^0 +SlJr). 



(5) 



The upper and lower states of the system are denoted as 
|-|-) and |— ) respectively. Its transition frequency is ujq. 
<Jz is the z-Pauli matrix in the basis {|-|-) , |— )}. Then, 
the contribution ThuqUz/'^ defines the bare system Hamil- 
tonian. The constants {Sujr} define the spectral shifts 
associated to each bath state. The second contribution 
in Eq. ^ introduces the interaction between the system 
and the external laser excitation, whose frequency is ujl- 



The operators cr^ — |+) (— | and a — |— ) (+| are the rais- 
ing and lowering operators acting on system eigenstates. 
The Rabi frequencies {flij} measure the strength of the 
system-laser coupling for each configurational bath state. 
The rest of the system operators appearing in Eq. ([Sj are 
defined by 



D = crV/2, 



J['] 



(6) 



while {•••}+ denotes an anticonmutation operation. 
Then, the contribution proportional to the constant jr 
defines the natural decay of the system associated to each 
reservoir state. 

The second line in Eq. ^ introduces a coupling (with 
rates rjjiipi) between all the states {pfl;(t)}, representing 
the fluctuations (transitions) between the configurational 
states of the environment. Depending on the definition 
of the system operator A different cases are recovered. 
When A = I, where / is the identity operator, the tran- 
sitions between the configurational states do not depend 
on the system state. Hence, the probabilities ([2]) are gov- 
erned by a classical master equation whose structure fol- 
lows straightforwardly from Eq. ([3]) . This case allows us 
to describe situations such as spectral diffusion processes, 
conformational environment fluctuations that affect the 
natural decay of the system, as well as single fluorophore 
systems diffusing in a solution [SH . When Aj^ I, the con- 
figurational fluctuations are statistically entangled with 
the state of the system. Depending on the structure of A 
different kind of situations can be described such as for 
example light assisted process, where the fluctuations of 
the bath depend on the external laser field intensity. 



Vectorial representation 

In order to establish a general formulation of the 
quantum-jump approach, we introduce a vectorial no- 
tation that allow to simplifying the presentation and cal- 
culations. To the configurational bath states we asso- 
ciate a vectorial space, defined by a basis {1^)}^;™, with 
{R\R') = 5rr', each vector \R) being related to a differ- 
ent configurational bath state [S^. The set of auxiliary 
states {p_R,(i)} allows us to define the vectors 

\p,) ^ 5]pH(i)ii?), \Pt) ^Y^^vswrnR)- (7) 

R R 

These two objects encode both the system dynamic and 
the evolution of the configurational bath states. In fact. 



ps{t) = (l|pt). 



PR{t) ^ {R\Pt) , 



(8) 



where we have defined the i?-vector (1| = '^r{R\- These 
identities follows straightforwardly from Eqs. ([1]) and ([2]) 
respectively. The normalization of the system state can 
be written as Tr5[(f |pt)] = 1, while the normalization of 
the configurational populations read (l|-Pt) = f . 



With the vectorial notation, the Lindblad rate equa- 
tion (J3|) can be rewritten as 



d\pt) 
dt 



^C\pt). 



(9) 



The structure of the matrix of system superoperators C 
follows from (J3]). From now on, with the hat symbol 
we denote vectors in the i?— space whose components are 
superoperators acting on the system Hilbert space. 



III. QUANTUM-JUMP APPROACH 

Our goal is to characterize the photon emission process 
associated to the fluorescent system. Of special interest is 
to determine how the environment fluctuations broke the 
renewal property in successive photon emissions. This 
property, for example, can be easily determine from a 
single experimental realization by measuring the succes- 
sive time intervals, {tj = ti — ii-i}, between consecutive 
photon recording events (happening at times U and ti_i). 
Then, one can define the waiting time distribution 



^W(T)^(5(r-r,)) 



real ' 



(10) 



where (• • • )j.^^j denotes a time average along a single real- 
ization [{f{n))^^^, = \imt^^il/t)J^dt'f{n{t'))]. Conse- 

quently, Woo (ti) defines the stationary probability den- 
sity of the intervals {ri}. It satisfies the normalization 

J^ drw oo {t) = 1. Similarly, one can define the (station- 

(2) 

ary) probability distribution Woo {t2 , ti ) for two consec- 
utive intervals (r^ and Ti+i), i.e.. 



,(2) 



(7-2, Tl) = ('5(t2 - Ti+l)5{Tl - Ti)) 



,(2) 



real ' 



(11) 



It fulfills /p dT2 J^ dTiWoo{T2,Ti) = 1, and the con- 
sistency relations J^ dT2Woo{T2,Ti) — Woo (ti), and 

/p°° dri w^^' (r2 , Ti ) = w^cx! (t2 ) . 

By knowing both probability distributions, one can 
quantify how much the photon emission process departs 
from a renewal one. The departure from zero of the di- 
mensionless parameter 



A(r2,Ti) 



(2),' \ 

_ Woo(T2,Ti) 



„(1) 



ICoo (T2)wio (ti) 



(1)/ 



(12) 



measures the strength of the non-renewal effects induced 
by the bath fluctuations. In fact, in absence of fluctu- 
ations the system dynamics becomes Markovian obey- 
ing the relation wL (t2,ti) = wL (r2)wL (n), imply- 
ing A(r2,Ti) = 0. Nevertheless, we remark that non- 
Markovian system dynamic may also lead to renewal 
emission process |33M35| . 

The possibility of flnding analytical expressions for 
wio (ti) and wio {t2,ti) are one of the central results 
of this contribution. We solve this task by extending the 
quantum-jump approach on the basis of Eq. ^ [Eq. ^]. 



A. Measurement operators 



Here, we have introduced the non-unitary propagator 



The quantum-jurnp approach rehes on a quantum mea- 
surement theory [l|-|3| . Here, the definition of a measure- 
ment operation must to include both the system and the 
configurational bath states. If \p) is the state previous to 
a measurement, the state AA^\p) after measurement is 



MM 



Mp) 



Trs[(l|Jjp)] 



(13) 



The vectorial superoperator J^^ define the unnormalized 
transformation of \p) due to the measurement action. 

When not any measurement is performed over the con- 
figurational space, one must to consider only one super- 
operator A^^ associated to the photon detection events, 
p — > photon- detector [see Eqs. (gU and (155)) ]. Neverthe- 
less, we will also consider the existence of extra measure- 
ment channels that may provide information about the 
configurational states of the reservoir [see Eqs. (IA2p and 
(IA10[) ]. Hence, Eq. (jH]) is discomposed as 



d\pt) 



dt 



= (23 + ^J^)|pO: 



(14) 



where T) 



L 



^ J'^. In the Markovian case, i.e., 
when the configurational space is one-dimensional, the 
(unique) superoperator J^ is related to the wave vector 
collapse after a photon recording event, while T) defines 
the conditional dynamic between consecutive photon- 
detections [il-Q- Here, the formalism must also to take 
into account the fluctuations of the environment, i.e., 
the vectorial nature of \pt) and the existence of different 
channels (labeled by p) that may also provide informa- 
tion about the transitions between the bath states. 

Equation p4t provides us the basis for characterizing 
the recording process. The following formulation is gen- 
eral, being independent of both the specific structure of 
Eq. ([3]) and the definition of the measurement channels 
{M.fj}. Specific examples are worked out in Section IV 
and Appendix A. 



f{t,T)\p)^ 



^■D(t- 



^^|P) 



Tr5[(l|e*(*-)|p)]' 



(17) 



the function 



Po[i,T;|p)]^Tr5[(l|e*(*--)|p)], 



(18) 



and the scalar contribution 



r^[|p)]^Trs[(l|J'Jp)]. 



(19) 



By associating the propagator T{t,T) with the (vec- 
torial) conditional system dynamic between consecutive 
recording events (photon-detections or/and configura- 
tional transitions), the first line of Eq. ([T6| can be in- 
terpreted as the contribution of all measurement realiza- 
tions where not any detection event happens up to time t. 
Consistently, the weight Po[t,T; \p)] must to interpreted 
as the corresponding (survival) probability for not hav- 
ing any transition in the interval (r, t) , given that the last 
one happened at time t, where system state is \p). 

The second line (integral term) of Eq. ([T6l) can be read 
as the contribution of all realizations where a measure- 
ment event happens at time r [represented by the ac- 
tion oi M^ on \pr)] and not any detection happen up 
to time i, which justifies the presence of T(t, r) and 
the survival probability Po[t,T;Mf_i\pT)]- Consistently, 
r^[|pr)]rfT must to define the probability of having an 
event in the ^-detector in the time interval (t, t + dr) . 



By expressing Eq. pB]) as a sum over all possible mea- 
surement outcomes, the previous statistical interpreta- 
tion can explicitly be demonstrated. By writing 



B. Statistic of the detection events 

The statistics of the successive recording events can be 
obtained after writing the system dynamics as an inte- 
gral over all possible measurement paths. The evolution 
Eq. P^ can formally be integrated as 



\pt) = e^* Ipo 



)+E 



.^(t-^: 



J^ \pr) dT, (15) 



which can straightforwardly be rewritten in terms of the 
measurement operators {Mfj,} as 

\pt) = Po[i,0;|po)]r(i,0)|po) (16) 

+y2[ Po[t,T;MM)]fit,T)MM)r^[\Pr)]dT. 
„ ■'0 



|p,) = ^(i)|po) = £^(")(t)|po), 



(20) 



n=0 



with ^(o)(t) ^ pp[i 0, \po)]f{t,0), from Eq. dUD we get 



^(")(i) = V / dt„ ■■■ f\ti Pr,[t, {t,}^, {//.iri 

xfit, t„)M^„ ■ ■ ■ fit2,ti)M^,,f{ti,0). (21) 



The weight P„[t, {t,}'/, {piji] is defined by 



Pn[t,{t^}'lA^^^}l] ^ Poit-.tn, M f,„ \ptj] w^,Jtn,tn-i; M f,^_,\pt„_,)] ■■■ w^^[t2,h; M ^,,\pt,)] Wf,,[ti,0; \po)]. (22) 



r 



The intermediate states read 



\Pu^^)^T{U+l,U)M^,^Pu), (23) 

with |ptj = f (ti,0)|po), while 

«;^[i,r;|p)]^Trs[(l|J,,e*(*-^V)]- (24) 

Clearly, Q^^^t) [Eq. (PT|) ] can be associated to all trajec- 
tories where n-detection events happen up to time t, each 
one at times {ii}*^" in the {pi}i detectors. The interme- 
diate evolution between detection events [A^/^J is given 
hyf{U,U.i). Consistently, P4t,{uy{,{^l,}^] [Eq. dH])] 
defines the probability density of each trajectory. Thus, 
Wp,i [ti,ti-i;A4fi,-_-^ \pti-i)]dti can be read as the probabil- 
ity of having a detection event in the /i^-detector in the 
interval {ti,ti + dti) given that the last detection event 
happened at time ti-i in the /ii_i-detector, not happen- 
ing any event inside the interval (ii_i,ii). 

By using the normalization of the vectorial state, 
(d/dt)Trs[(l \pt)] = 0, from Eq. ^ it follows the rela- 
tion Tr5[(l|P|.)] = -E^TrsKllJ-J.)]. Then, Eq. m 
can alternatively be written as 



,, Jo 



\p)]dT. (25) 



With this relation, we notice that Eq. (1221) has the 
same structure than a renewal process, i.e., there ex- 
ist a probability distribution (waiting time distribution, 
Wp,[t,T;\p)]) that define the statistic of the time inter- 
val between consecutive detection events. Nevertheless, 
here the waiting time distribution has a functional de- 
pendence on the system state posterior to a detection 
event [A^^jp)], which broke the renewal property. 
By writing the states \pti+i) [Eq. ([M])] as 



\PU+i 



s*(*'+^-*')j,JptJ 



(26) 



Trs[(l|e^(*.+i-*')J'^JptJ] 
the n-joint probability density (j22l) can be rewritten as 



Pn[t,{t^y^A^^^}l] = Trs [(1 [e^^*-*") J'^„ • • • J-^, e^^*^"*^ ) 

(27) 



X J-^.e^*^ \po)] 



This expression recovers the result of Ref. [331 . Notice 
that its structure is similar to that obtained in the con- 
text of a photon measurement theory [l|-[3|- Nevertheless, 
here the underlying trajectories are vectorial and depend 
on the extra parameters pi, i — 1 ■ ■ ■ n. 

The probabilities P„(t) of having n-detection events 
up to time t can be obtained by integrating the joint 



probabilities densities P„[i, {ii}" , {pi}i] over all possible 
detection paths 



Pn{t) 



Mil ■■■Ml 



dtr, 



dhPn[t,{uy,{p,y]. 



(28) 

In the context of SMS, objects of this kind are usually 
characterized through a generating function approach 
based on a stochastic Bloch equation |26l - l3G |. Then, 
while previous approaches are able to get these objects, 
the present treatment also allows us to get the underlying 
joint statistic defined by Eq. (P^ . 



C. Stationary waiting time distributions 

The joint probability density Eq. (|22|) is one of the cen- 
tral results of this section. It completely characterizes 
the statistic of the recording events. It can experimen- 
tally be determine from an ensemble average over mea- 
surement realizations having n-detection events in the 
interval (0,i). Nevertheless, the stationary waiting time 
distributions Eqs. (jlOl) and pT|) are defined by a time av- 
erage along a single measurement realization. For ergodic 
environment fluctuations, objects of this nature can be 
studied by describing the measurement process after hap- 
pening an infinite number of recording events and that an 
infinite time elapsed since the initial condition, \po) . In 
Appendix B, we show that in that limit Eq. (j22[) remains 
valid under the replacement 

\po)^M\Poo), (29) 

where |/0oo) corresponds to the stationary state 



= lim \pt)- 



(30) 



It comes forth because a time averaging procedure can 
only provides information about stationary observables. 
The measurement operator Ai is defined by 



M\p) 



J\P) 



Tr5[(l|J|p)] 



^^E^A^' 



(31) 



and takes into account the happening of an arbitrary 
measurement event in the long time regime. With these 
definitions, from Eq. ([^^ we introduce the first station- 
ary waiting time distribution 



.,(1) 



(t,p) = w^[t,0;7W|poo)] 

= Trs[(l|J'^e*^>l|poo)] 



(32) 



as well as the second order stationary waiting time dis- 
tribution 

W^^Kt2,H2;Ti,Hi) = W^^[t2 +Ti,Ti;Mf,,\pr,)] (33) 

xWf,,[Ti,0;M\p^)] 

Higher objects, Wc» [{T,;}",{/ii}"], can be written in a 
similar way. They define, in the stationary regime, the 
probability density of the time intervals JTi}" between 
successive recording events happening in the {/Xi}" detec- 
tors. When the measurement process only involves a pho- 
ton detector apparatus, ji —^photon-detector, Eqs. (l32)) 
and (1551) allow to get analytical expressions for the dis- 
tributions ([TU]) and ([TT|) respectively (see Section IV). 



D. Stochastic density matrix evolution 

From the previous analysis, we obtained the recording 
event statistics associated to the density matrix evolution 
Eq. (J14p . The quantum-jump approach also allows build- 
ing up the underlying stochastic dynamics that reproduce 
that statistic. The key ingredient is the definition of a 
stochastic process developing in the system Hilbert space 
and whose realizations can be mapped with the realiza- 
tions of the measurement apparatus signals. The average 
over realizations must to recover the system density ma- 
trix evolution. Then, in the present context we search 
for the definition of a stochastic vector \pf), such that 
\pf) = \Pt), where \pt) is defined by the evolution ([T4l). 
From now on, the overbar denotes (ensemble) averaging 
over realizations. 

Based on the path integral solution obtained previously 
[Eq. (|20l)]. the stochastic evolution can be written as a 
piecewise deterministic processes Q 



,t,dNt 



j^\pf) = [T>-Trsil\V\pT)]\pf) + J2i^,~l)\pf) ^^ ■ 

M 

(34) 
Here, the deterministic non-linear term [first contribu- 
tion on the r.h.s.] corresponds to the conditional evolu- 
tion in the intervals between consecutive measurements 
events, i.e., the dynamics defined by Eq. pT)) . On the 
other hand, the second term introduces the disruptive 
changes in the vectorial state after a measurement event, 
i.e., \pf) — 7> A^plpf ). Consistently, the noisy terms are 
defined by dN^/dt = J2k^i^ ~ ^k)^ where t^ are the 
times where a measurement event happens in the p- 
detector. By denoting with Nj^ the number of detections 
events up to time t, it follows the alternative definition 
dN^ — {N^,^^ — Nj^), i.e., dNj^ are the increments of the 
(Poisson) process Nj^ [3] . In agreement with the previous 
analysis, their average must to recover Eq. ([TO| . i.e.. 



By using the property dNidN]^ — S^^idNi , which im- 
plies that a simultaneous detection in two different mea- 
surement apparatus is never observed and that (dNj^)'' = 
dNi ^ in Appendix C we show that Eq. (fT4|) is recovered 
after averaging Eq. ([M)) over realizations. 

The realizations associated to Eq. ([M)) can be eas- 
ily determine after providing a recipe for calculating the 
random times where the detection events happen. Their 
numerical calculation relies on evaluating the statisti- 
cal objects introduced in Eq. (IT51) along each trajectory. 
Given that the system is in the state |pf ), the quan- 
tity F ti[\pf)]dt [Eq. (Uni)] gives the probability of having 
an event in the /i-detector in the time interval (i, t -\- dt). 
This quantity defines an infinitesimal time step algorithm 
(see Appendix D). In a similar way, P[)[t^t'\M.fj\pf;)] 
[Eq. (dHI)] define the survival probability for the next de- 
tection event (at time t) given that a /x-detection event 
happened at time t' . This object allows to defining a fi- 
nite time step algorithm (see Appendix D). 

Independently of the method (algorithm) used to de- 
termine the times of the recording events (transitions), 
given that at time t a measurement happens, \pf') — >■ 
M.ij\pf), each transformation TM^ [Eq. (fT3|) ] must be 
chosen with probability 



V(i) 



FMt)] 



Tr5[(l|JMlpf)] 



Y.,'rA\pf)] E,' Tr5[(i| j^Hpf)]^ 



(36) 



which satisfy E/it/i(i) = 1- This rule corresponds to a 
selective measurement of the set of ^-observables Q . Be- 
tween successive recording events, the evolution of \pf) 
is deterministic and defined by Eq. (IT71) . 
Through the relations 



p^S*(i)^(l|pf), 



Wf) 



Tr5[|pf)] 



(37) 



the vectorial state \pf) provide a stochastic representa- 
tion of both the system density matrix [Eq. ([Ij], pg (t) = 
ps{i)i and the occupation of the configurational bath 
states [Eq. ©], {R\Pf) = PR{t). In contrast with the 
standard quantum-jump approach, in general it is not 
possible to get a simple dynamical evolution for p^g{t) 
[or to (i?|P|')]. In fact, here the formalism relies on the 
vectorial nature of \pf) [however see also Appendix A]. 



E. Non-renewal recording realizations 

The trajectories associated to Eq. (p4)) allow us to es- 
tablishing a simple scheme for understanding the non- 
renewal nature of the recording process. In fact, its un- 
derlying structure is similar to that of a renewal one. 
Given that the last event happened at time t' in the p- 
detector, the random time t for the next event is defined 
by a waiting time distribution Wst{t,t' , p), which read 



Wst{t,t',p) 



dt 



Po[t,t';M,.\pn, 



(38a) 



dN,^ = TTs[{l\J^\pt)]dt = r^[\pt)]dt. 



(35) 



-TTs[il\Ve*^'-''^M^\pp)]. (38b) 



By using the relation Trs[(l|P|.)] = - E^ Trs[(l|J-^|.)], 
it follows 

wstit,t',fi) = Trsmje*('-''^M^\pn, (39) 

where J = J2ij,Jf^ [Eq- (EH)]- At time t, \pf) is up- 
dated with the conditional evolution Eq. ([T7|) and the 
new recording event is selected with the probabilities 
([55)1 . The next events follow from the same rule (see 
Appendix D). The average over realizations recover the 
statistics defined by Eq. (1^ . 

The departure of the recording realizations with re- 
spect to a renewal process comes from the dependence 
of Wst{t,t' , p) on Ipj'). Only if M^\p^^) is independent 
of |p(*) one get a renewal recording process. Never- 
theless, in general this does not happen, implying that 
Wstit,t' , fi) randomly change between successive events. 
Then, in contrast with a renewal process, here the suc- 
cessive events are defined by a stochastic waiting time 
distribution that parametrically depends on the vectorial 
state Ipj*). Finally, we notice that Wst{t,t' , p) can consis- 
tently be written as Wst{t,t',p) = I];i "'/i[^'*';-^A'l/°t*)]' 
where w^[i,i';p)] is defined by Eq. pi)) . 



IV. PHOTON EMISSION MEASUREMENTS 

In the previous sections, we developed a general theory 
that allows to characterizing the measurement processes 
associated to a broad class of physical situation arising in 
SMS. The theory depends on which kind of measurement 
process is performed over both the system and the config- 
urational states. In this section, we analyze the situation 
where there exists only one measurement process defined 
by a photon detector apparatus coupled to the scattered 
electromagnetic field. This is the standard situation in 
SMS, where any direct information about the configura- 
tional space is unavailable. Then, the parameter fi only 
includes one term corresponding to the photon detector. 
Furthermore, our formalism is able to describe different 
kind of environmental fluctuations. First, we analyze the 
case of self-fluctuating environments, i.e., when the tran- 
sitions between the configurational states do not depend 
on the state of the system. As a second leading case, we 
analyze environmental fluctuations that depend on the 
intensity of the laser excitation. 



given by 

^^fl.(i) ^-Y. ^R'RPRit) + Y. ^RR'Pr'W- (41) 

R' R' 

Hence, the stochastic dynamics between the configura- 
tional states is governed by a classical master equation 
that does not depend on the state of the system. This 
case allow to describe processes such as spectral fluctua- 
tions, life time fluctuations, and molecules diffusing in a 
solution. 



Photon measurement operator 



The measurement operator, Eq. (J13p . must to take into 
account all contributions that, independently of the R- 
state of the reservoir, lead to a photon emission. Then, 
from Eq. (gOl), we write {p, -> ph) [\p) = Y^j^ \R)pr] 



-^Phlp) 



Jphip) 



^ J2rJr\R) ^pr(^'^ 

Trs[(l|iph|p)] Efl' lR.'TTs[<T^apR, 



(42) 



Notice that each contribution in the sum corresponds to 
the standard definition arising in Markovian fluorescent 
systems [l-Q, i.e., A^phP = crpcrVTrsifVp]. The vecto- 
rial superoperator 2? [Eq. (|14p ] here is defined from 



V = C-J^ 



phi 



(43) 



where C follows from Eqs. ([9]) and (|40p. while ijph from 
Eq. (|42|) . Alternatively, V can be explicitly defined 
through the non-unitary evolution generated by it 



-(i?|pn = {R\V\pY), 



(44) 



where the index u say us that the auxiliary vector |p") 
is not normalized to one. In fact, its norm is related to 
the survival probability Eq. (fT8|) . By denoting p\{t) = 
{R\p't), we get 



dt 



'-[Hn,pl{t)]-^R{D,pl{t)}, 



E' 

R' 



'A'A/AW + E-^fi^' '"«'(*)• (45) 



Having the definition of the superoperators A^ph and T) 
we can apply the theory developed in previous section. 



A. Self-fluctuating environments 

This case is covered by Eq. ([3]) by taking A ~ I, 
dpR{t) -i, 



dt 



^ [HR,pR{t)]--iR{{D,pR{t)}+ - J[pR{t)]) 
J2'l'R'RPR^*^ + J2'^RR'PR'^^^- (40) 



For notational consistency we take riRR/ -^ (j)RR' W^ ■ 
From Eq. (HUl) . the evolution of populations Eq. ^ is 



2. Stochastic dynamics 

The dynamic of the stochastic state \pf) follows from 
Eq. (|M)) . The disruptive transformation associated to 
a photon detection event, |pf ) — > A^ph|pf ), from the 
expression (|^ . and by using that a = |— ) (+| , a^ = 
1+) (— I , can explicitly be written as 



\pf)^M^,^pf) = \-){-\Ypt{t)\R). 



(46) 
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FIG. 1: Stochastic realizations of a fluorophore system de- 
fined by the evolution Eq. (|40p . The configurational space 
is two-dimensional, R = A,B. The parameters are ^r = fi, 
SiOR = 0, 7A/n = 1, 7B/n = 10, 0AB/n = 0.003, 0flA/n = 
0.009. The laser is in resonance with the system, i.e., ljl ~ too- 
(a) Realization of the of the upper population of the system 
(+1 p'si^) 1+) ■ (b)-(c) Realization of the configurational popu- 
lation of the bath {R\Pf^), for R = A. (d) Intensity realization. 
The values Ir are defined by Eq. (f50)) . 



Here, the weights {p%{t)} satisfy the normalization 
J^RPRi^) = Ij ^i^d are defined as 



pm 



7R{ + \pUt)\+) 
ER'lR'{ + \P%it)\ 



(47) 



where the notation pfj(t) = {R\pf) was used. From 
Eqs. ([37| and (|46l). it is simple to get 



(48) 



pf{t)^{l\M^^\pf) = \-){-\, 



and that 



(i?|Pf ) ^ Tr5[(i?|^ph|pf )] = Ptit). (49) 

Eq. (|48)) shows that in fact, after a photon detection event 
the system collapses to its lower state |— ) . On the other 
hand, Eq. (|49|) say us that p%{t) is the value of the con- 
figurational populations after a photon recording event. 



Between the detection events the stochastic dynamics 
is defined the conditional evolution defined by the super- 
operator V [Eqs. (|43l) and (|45|) ]. On the other hand, as 
there exist only one measurement apparatus, the weights 
{t^{t)}, Eq. dMl), here reduce to tph(i) = 1. 

In the next figures, we consider a fluorophore sys- 
tem coupled to an environment characterized by a two- 
dimensional configurational space, R = A, B, which only 
affect the decay rates {7^} of the system, i.e., the Rabi 
frequencies [Eq. (JU] do not depend on the configurational 
states, fin = H., and the spectral shifts [Eq. ([5])] are null, 
dujR = 0. Furthermore, the laser is in resonance with the 
system, i.e., ljl = i^q- 

In Fig. 1(a) we show a realization of the upper pop- 
ulation of the system {+\pfit)\+) [Eq. ^]. The re- 
alizations were determined by using the finite time step 
algorithm defined in Appendix D. Each collapse of the 
upper population to zero is related to a photon emission 
[see Eq. (p|]. 

In Fig. 1(b), we show the realization of the configu- 
rational population {R\Pf) ior R = A [Eq. dST])]. As 
the configurational space is two-dimensional, {A\P^^) + 



{B\P! 



1. We remark that these realizations are as- 



sociated to a measurement process that only gives in- 
formation about the photon emission events. Not any 
information is provided about the configurational states 
of the bath. Therefore, the realizations of {R\P^'^) are 
the best estimation 40] about the configurational state 
of the reservoir that can be obtained by knowing the mas- 
ter equation ()40p and a given realization of the photon 
detector apparatus. 

In Fig 1(c) we plot {R\P^^) (for R — A) over a larger 
time interval. For the chosen parameter values, the con- 
figurational populations develop a quasi-dichotomic be- 
havior. When {R\P^^) « 1, we can affirm that is highly 
probable that the bath is in the configurational state \R). 

In Fig. 1(d), we plot the scattered intensity, which is 
defined by I{t) = [n{t -\- St) — n{t)]/St, where n{t) is the 
number of photon recording events up to time t and 6t 
is an adequate time flag averaging. Its fluctuations are 
highly correlated with the values of {R\Pl*^). In fact, the 
intensity fluctuates around two well-defined values Ir, 
which are defined by the intensity of a Markovian flu- 
orescent system [j-jj] characterized by the jgarameters 
corresponding to each configurational state [31j . i.e.. 



/r = 



IB^^R 



ll- 



2n% 



45|' 



(50) 



where Sr = ujl — ujr. The dichotomic behavior arises be- 
cause the system is able to emit a large number of pho- 
tons previously to the occurrence of a transition between 
the configurational bath states, i.e., X]fl' (f'R'R ^ Ir- 

In Fig. 2(a) we plot the upper population (-|-| ps{t) |+) 
that follows from Eq. ([^O)). as well as an average over 
(«103) realizations of (+1/3^(01+) [see Fig. 1(a)]. In 
Fig. 2(b), we plot the analytical solution of the configu- 
rational populations {R\Pt) defined by Eq. (HJ), as well as 
an average over realization of {R\P^^) [see Fig. 1(b) and 
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FIG. 2: Evolution of the upper population (+|ps(f) |+) (a) 
and the configurational populations [R\Pt) (b), that follows 
from Eqs. (|40p and (|4ip respectively. The parameters are the 
same than in Figure 1. The initial conditions are ps(0) = 
|-) (-j and (AjPo) = (B\Pq) = 1/2. The noisy curves follow 
from an average over the realizations shown in Fig. 1. 



(c)] . In both cases the ensemble averages recover the dy- 
namics dictated by the corresponding master equations, 
showing the consistency of the developed approach. 



3. Photon emission process 

The recording events are characterized by the stochas- 
tic waiting time distribution Eq. p9p . Then, we write 

w.t{t,t') = Trs[(l|ie*(*-*')^ph|p?^)]. (51) 

Notice that here J' — J7ph. The function Wst{t,t') define 
the statistic of the time intervals between consecutive 
photon emissions. From Eqs. (|42l) and (H51) it follows 



Wst 



{t,t') = j2^R' 



-RR/ 1 



\PR'it')]\+). (52) 



RR' 



-T} ' r 

p;(t')=i.o 

p;'(f)=0.75 

p;'(f)=o.5 

PR(f)=0.25 
PrCH^O.O 




This expression allows us to get an analytical expression 
for Wst(t, t') [not provided due to its extension] that para- 
metrically depends on the set {p%{t')}, Eq. (H7|) . 
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FIG. 3: (a) Plot of the waiting time distribution Wst{t,t') 
[Eq. ((52])] for different values oi p%{t'). From top to bottom 
we take {R = A) pt{t') = 1, 0.75, 0.5, 0.25 and 0. (b) Sta- 
tionary waiting time distribution [Eq. (1531) ]. The stationary 
weights [Eq. ^] readpX = 0.527 andpg= = 0.473. The noisy 
curve corresponds to a numerical distribution determine from 
the intervals between consecutive photon emissions [Eq. (I10|) ] 
along a single realization (like that shown in Fig. 1). In both 
plots, the parameters are the same than in Fig. 1. 



As the set of weights {p^, {t')} correspond to the config- 
urational populations after a photon recording event [see 
Eq. (I49p ]. the waiting time distribution change between 
consecutive photon emissions. The successive (stochas- 
tic) values of {p%{t')} can be read from the realization 
of (i?|Pjf') shown in Fig. 1(b). Notice that in each event, 
defined by the collapses (+|/3s(*)l+) ^ [Fig- 1(a)], 
{RlPf*") suffer an abrupt change in its slope. 

In Fig. 3(a), we plot Wst{t,t') (as a function of t — t') 
for different values of the parameters p^, (i'). As the con- 
figurational space is two-dimensional, the two parame- 
ters satisfy the normalization p^(i') + p%{t') = 1. We 
notice that Wst{t, t') has a strong dependence on the val- 
ues of the configurational populations {p%{t')}, which 
in turn say us that the photon emission process strongly 
departs from a renewal one. For Markovian fluorescent 
systems, the set {p%it')} reduce to only one parame- 
ter with value equal to one (the configurational space 
is one-dimensional). Therefore, Wst(t,t') is the same ob- 
ject along a measurement trajectory, recovering a renewal 



10 




ioI\t2,ti)I^^ 



Qx 



1.0- 



0.5- 



<' 0.0- 



-0.5- 



"(b) 


1 


' 1 


' 1 ' 




/^T=15 


■ 


Qt,=5 


■ 


- \_^ 








:^ 






^\=^ 


■ 


Qt=3 







1 
5 


1 1 

10 
Qx 


' 1 1 

15 


2 



FIG. 4: (a) Stationary two-time waiting time distribution 
Wio (r2,ri) [Eq. (|56p ] for difTerent values of ri. From top 
to bottom we take Q.t\ = 3, 1, 5, 0.3, and 15. The noisy 
curves correspond to the time average Eq. (Illf) . (b) Parame- 
ter A(t2,ti) [Eq. dni)] determine from Eqs. ^^ and ^^, for 
different values of n. From top to bottom Q.t\ = 15, 0.3, 5, 
1, and 3. The parameters are the same than in Fig. 1. 



process. 
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FIG. 5: Contour plot of w\j {t2,t\)/Q? [upper panel] and 
A(r2,ri) [lower panel] corresponding to the plots of Fig. 4. 



4. Stationary waiting time distributions 

By measuring the time intervals between successive 
photon emissions along a given trajectory one can deter- 
mine the stationary waiting time distribution Eq. (|10p . 
The analytical expression for this probability distribution 
can be read from Eq. (|32p . We get 

w'^^\r) ^Y.^R (+1 e%[\~) (-[p^,] 1+) , (53) 

RR' 

where the weights p'^ are defined from the relation 

K^\Poo) = \-){-\Y.p'^\Rl (54) 



delivering the expression 

_ ir{Mpr 



Pr 



Er'1R'{+\p^'\+)' 



(55) 



Here, p^ = limt^oo PRit) = {R\poo) [Eq. iOD]. By com- 
paring Eq. ((53)) with Eq. ([52|) . we realize that wL (t) fol- 
lows from Wstit, t') after the replacements p%{t') — >■ p^, 
and {t — t') -^ T. 

In Fig. 3(b) we plot the analytical expression for 

Woo {t) that follows from Eq. (|53p . Furthermore, we show 
a numerical distribution determine from the time aver- 
age Eq. (fTU]) . The theoretical distribution correctly fit the 
numerical result. Consistently, we also checked that the 
time average of {p%{t')} along a single trajectory recover 
the weights {p^}, Eq. ((55|) . 

The analytical expression for the second waiting time 
Eq. pT]) can be obtained from Eq. (I55|) . delivering 



.(2)/ 



WX^'{t2,Ti) = 2_^lR 
RR' 



-RR' I 



kff(n)] 



(56) 
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where the functions (^jj(ti) read 



iPr{ti) = 7_r, V(+|e 



R' 



T>Ti 

RR' 



I 71°° 1 

\Pr'\ 



(57) 



In Fig. 4(a) we plot the analytical expression for 
wL (t2,ti) that follows from Eq. ([56)) for different val- 
ues of Ti. Furthermore, we show the numerical result 
that follows by determining the probability distribution 
of two consecutive time intervals between successive pho- 
ton emissions along a single trajectory, Eq. (TTTI) . The 
theoretical result correctly fit the numerical distribution. 

In Fig. 4(b) we plot the dimensionless parameter 
A(T2,ri) [Eq. ((H])] determine from Eqs. 1^ and (I5S1) . 
for different values of ti . For almost all values of T2 and 
Ti, A(r2,Ti) departs appreciably from zero, indicating 
the departure of the photon emission process from a re- 
newal one. We also note that there exist special values 
of the consecutive time intervals where A(t2,ti) ~ 0. 
From the definition Eq. (IT^ . we deduce that when 
A(t2,ti) > 0, the frequency of the successive intervals 
Ti and T2 is greater than in the renewal case. The situa- 
tion A(t2,ti) < 0, admits the inverse interpretation. For 

clarifying the structure of both Woo (t2, ti) and A(t2, ti), 
in Fig. 5 we show their contour plots. As can be deduced 
from Fig. 4 and 5, A(t2,ti) reach it maximal values for 
higher values of both T2 and n . 



5. Slow and fast environment fluctuations 

The expressions for the stochastic waiting distribution 
Wstit,t') [Eq. dSg)], and the first [w^^\t), Eq. ^] and 
second [woJ (t2, ti), Eq. ([55]) ] stationary waiting time dis- 
tributions allow to characterize the photon emission pro- 
cess as well as its departure with respect to a renewal one. 
Here, we provide simple analytical expressions for these 
objects in the limit of both fast and slow environment 
fluctuations. 

The characteristic time of the bath fluctuations are 
measured by the rates {(J)r'r}, Eqs. (|40|) and (|4T|) . On the 
other hand, the average time between photon emissions is 
measured by the inverse of the intensities {Ir}, Eq. ([5n| . 

When the bath fluctuations are much slower than the 
average time between photon emissions, {(f>R'R} <C {/_r}, 
it is valid to approximate the conditional evolution de- 
fined by the superoperator V [Eqs. (j331) and (01])] as 

4r^'\H HPR'it')] « SRR,elt''\\-) {-\PR{t')]- 

(58) 
This approximation corresponds to disregarding the non- 
diagonal contributions between photon recording events. 
By inserting this condition in Eq. (|5^ we get 

w.,itX) ^ T.^R{+\etV\\-) {-\pW)]\+) , 

R 

^ J2y^Rit-t')pt{t'). (59) 



Then, Wstit, t') can be written as a linear combination of 
the waiting time distributions {ii;^(i)}, each one being 
defined by the expression 



WR(t)^^R{+\e^^j,[\-){-\]\+) 



(60) 



This function correspond to the waiting time distribution 
associated to a Markovian fluorescent system [l|, Q with 
decay rate 77?, and whose Hamiltonian is given by Hr^ 
Eq. Q, i.e., its transition frequency is uir = (wq + 5ijJr), 
and its coupling to the external laser is measured by ^r. 
This result can straightforwardly be read from Eqs. (pi)) 
and (|151) under the replacement (pRR' — )■ 0. In the Laplace 
domain, t — > u, it can be written as [34| 



wr{u) = 



7/?72 



nyiR{u) 



u + jr/2 \u'^ + u^R + Q'j^hR{u) 
where the auxiliary function hR(u) is 

, , . [u + lR/2f 

{u + 7fl/2)2 + 6j^ 
and Sr ~ ujl — i^R- After Laplace inversion, we get 



, (61) 



(62) 



WR{t) = 



2^Rnj 



a 



exp(-7fli/2)[cosh(^+t) - cosh(C^i)], 

(63) 
where C| = bl - ml + ^l) ± Cfl.]'/V(2V2), with Cr = 

{b% + ml + ^l)]'- 167^^fl}'/'- When Sr = 0, the 
expression of Refs. |l|, Q is recovered. 

We have checked that Eq. (|59l) joint with Eq. (|63l) pro- 
vide an excellent approximation to the exact functions 
plotted in Fig. 3(a). On the other hand, Eq. (f58| is also 
useful for approximating the stationary waiting time dis- 
tributions. Eq. ([55)1 leads to 



W(r)c.^«;«(r)p^, 



while from Eq. ([5^ . we get 



w, 



(2; 



(T2,n) 



R 



wr{t2)wr{ti)p°j_ 



(64) 



(65) 



Furthermore, under the hypothesis of slow fluctuations, 
from Eq. (|i(7| we can approximate 7_r(+|p^|+) — 
IrP^. The constants Ir [Eq. ([50]) ] are the intensities 
associated to each configurational state R. On the other 
hand, P^ are the stationary values of the configurational 
populations Eq. (gT]), i.e., P]f = limt^oo -Pr(0- There- 
fore, from Eq. (|55|) we get the approximate expression 



Pr 



IrPr 



Er'Ir'P^ 



(66) 



Eqs. ([64l) . (|65]) . and (|66)) . also provide an excellent 
approximation to the exact analytical results plotted in 
Figs. 3, 4 and 5. Moreover, they have a clear physical 
meaning. In the slow limit each bath state establishes an 
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intensity regime defined by Eq. ((50)) . Hence, the statis- 
tic of tlie non-renewal plioton emission process follows 
from an average of the renewal statistic associated to 
each state [defined by wr{t)]. The weight of each contri- 
bution is p^. Consistently, this factors, which are the av- 
erage configurational populations after a detection event 
[Eq. ([5^ ]. are proportional to the intensities Ir and the 
stationary populations P^ related to each bath state. 

When the bath fluctuations are much faster than the 
average time between photon emissions, {4>b.'r} 3> {^_r}, 
the configurational populations reach their stationary 
values, P^ = liuit^oo PR{t), before happening many 
photon emissions. Hence, the fluorophore behaves as 
a Markovian fluorescent system with decay rate 7 = 
Y^rIrP^^ Riibi frequency Q. = Y^r^rPr^ and de- 
tuning 5 = J2r^rPr- The photon emission process 
becomes a renewal one, being defined by the waiting 
time distribution Eq. (p51) with {7/?, Q,fi, Sfi} — > {7, ft, S}. 
Near of this limit, for two-dimensional configurational 
spaces. Woo {t2, ti) develops small asymmetries on its ar- 
guments Woo (T2,ri) 7^ wio (ti,T2). In general, this prop- 
erty may arises in the intermediate regime between fast, 
Woo{t2,ti) ~ Woo (''■2)^00 (ti), and slow bath fluctua- 
tions, Eq. (1551) . 

The configurational fluctuations are frozen when 
Eq. (j40|) is defined with I^r'rI — 0. This case was par- 
tially addressed in Refs. [3J,|35|. Our present treatment 
provides a general description. Evidently, the configura- 
tional populations remain unaffected during all the evo- 
lution, \Pt) = |Po). The resuhs of Ref. [H [H follows 
from the approximation {P^^) ~ \Po): which is valid in a 
weak laser intensity regime and when the dynamics de- 
velops two different times scales induced by an infinite 
dimensional configurational space. 



B. Light Assisted environment fluctuations 

The general evolution Eq. ([3]) may also cover the case 
in which the statistical properties of the radiation pat- 
tern, as well as the environment fluctuations, depend on 
the external laser intensity [3l|, |3g|, i.e., light assisted 
processes. By taking A = a, and rjufj/ — ?> Jrr', we write 

^^ = Y^Hn,PRm-lRi{D,PRit)}+-J[pRim 

-J2iR'r{D, PRit)}+ +J2lRR'J[pR'{tM67) 



where D and J^ follows from Eq. ([6]). In this case, 
the evolution of the configurational populations [Eq. ©] 
strongly depend on the state of the system. In fact, here 
the configurational transitions may only occur when a 
photon emission happens. Thus, in general it is not pos- 
sible to write a simple equation for their evolution. Only 
when {^RRi} <C {7_r}, a classical rate equation similar to 
Eq. (gl]) can be derived (HHi. 



1. Photon measurement operator 

Here, the photon-measurement superoperator 
Mph\p)^ Jph|p)/Trs(l| Jphip), from Eq. §7}, reads 



■^Phlp) 



J2r \R){lR'^PR<^'' + J2r' 1RR"^PR'<^'^} 
Ej?"7i?."Trs[crtcrpfl"] 



where \p) = ^^ \R)pR^ and we have defined the rate 



(68) 



1R = 1R + X!fl/ ^R'R- 



(69) 



As in the previous case [Eq. (|42p ]. Eq. ((68|) take into 
account all possible configurational paths that lead to 
a photon emission. The conditional evolution defined 
by the operator T) = C — i/ph, expressed through the 
evolution of the unnormalized state |p") [Eq. (04])] reads 

^^ = ^[Hr, pim ^n{D, pl{t)}+. (70) 



Notice that in contrast with Eq. p5|) . here the conditional 
evolution is diagonal in the i?-space. 



2. Stochastic dynamics 

The structure of the stochastic dynamics of the vecto- 
rial state |pf ), Eq. ([34| . is similar to that of the previous 
case. When a photon detection event happens it implies 
the transformation 



\pf) ^ >Jph|pf 



\Y.ptmR), 



(71) 



where the weights satisfies 'YIirPr 
here they read 



1. From Eq. 



vtit) 



IR (+1 PRJt) 1+) + Y.R' IRR' ( + 1 PR'(t) 1+) 



j:R"lR"{+\P%'it) 



(72) 



From Eq. ([37l) it is simple to demonstrate that Eqs. 
and (H91) are also valid in this case. Therefore, in each 
photon recording event the system collapse to its ground 
state while p^it) define the posterior value of the config- 
urational populations. 

In the next figures we consider a fluorophore sys- 
tem whose evolution is defined by Eq. ([67| and a two- 
dimensional configurational space, R = A,B. The Rabi 
frequencies [Eq. ^] do not depend on the configurational 
states, J7fl = fl, and the spectral shifts [Eq. ([5])] are null, 
6uJu = 0. Therefore, the bath states only affect the decay 
rates {7^} of the system. The laser is in resonance with 
the system, i.e., wz, = uiq- 

In Fig. 6(a) we show a realization of the upper popu- 
lation of the system {+\ p^'^i^t) \+) [Eq. dST])]. The times 
of the photon emission events correspond to the collapse 
of the upper population to zero. Fig. 6(b) shows the re- 
alization of the configurational population of the bath 
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FIG. 6: Stochastic realizations of a fluorophore system de- 
fined by the evolution Eq. (|67|) . The configurational space is 
two-dimensional, R = A,B. The parameters are Qr = SI, 
Slor = 0, 7A/n = f.8, 7s/n = 0.15, 7As/n = 0.35, 
jba/^ = 0.2. The laser is in resonance with the system, 
i.e., ljl ~ ijJo- (a) Realization of the of the upper popula- 
tion of the system (+|ps(^) 1+) • (b)-(c) Realization of the 
configurational population of the bath (i?|P/'), for R = A. 
(d) Intensity realization. 7oo is defined by Eq. (|73[) . 



{R\Pt), ior R^ A [Eq. ^]. Due to the chosen param- 
eter values, at any time it is not possible to predict with 
total certainty [{R\Pl*^) w 1] the configurational state of 
the reservoir. This fact is evident from Fig. 6(c), where 
we plot {R\P^*^) (for R = A) over a larger time interval. 
In Fig. 6(d), we plot the scattered intensity. Consistently 
with the behavior of {R\P^*^), the intensity does not de- 
velop any dichotomic behavior. The intensity fluctuates 
around the value loo defined by (see Eq. (55) in Ref. |3l| ) 



R 



\Pr 



(73) 



where as before p^ = limt__j.oo Pi?(0- 

In Fig. 7(a) and (b) we plot the analytical solutions 
of the upper population {+\ps{t) \+) and the configu- 
rational populations {R\Pt) that follows from Eq. ((67|) 
[and Eq. ([2])]. The noisy curves correspond to an aver- 
age over realizations like those shown in Fig. 6. Notice 
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FIG. 7: Evolution of the upper population {+\ps{t) \+} (a) 
and the configurational populations {R\Pt) = Tts[{R \pt)] (b) 
that follow from Eq. (|67[l . The parameters are the same than 
in Fig. 6. The initial conditions are ps(0) — |— ) {— | and 
(AjPo) ~ {B\Po) = 1/2. The noisy curves follow from an 
average over the realizations shown in Fig. 6. 



that here, the behavior of the configurational population 
strongly depart from an exponential one, indicating that 
their underlying dynamics is highly non-Markovian. The 
physical origin of this characteristic is the dependence of 
the configurational transitions on the system state. 



3. Photon emission process 

In this case it is also possible to define a stochastic 
waiting time distribution that parametrically depends 
on the configurational populations after a photon detec- 
tion event, i.e., w,tit,t') = Trs[{l\Je*^*-*'^Mpi,\pf)l 
Eq. (151|). From Eqs. ^8^ and ((HI) we get 



Wst 



(t,i') = E^«' 



-RR' 1 



\PR'{t')] 



(74) 



RR' 



where p%{t') is given by Eq. (|72|l . As the vectorial su- 
peroperator T) is diagonal [see Eq. (jTOI ]. this expression 
can be written as 



Ustit,t')=J2wRit-t')p%{t'), 



(75) 
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FIG. 8: (a) Plot of the waiting time distribution Wst{t,t') 
[Eq. ((751)] for different values of p%{t'). We take the val- 
ues Pfl(t') = (full hue), 0.7 (dotted line), and 1 (dashed 
line), where R = A. (b) Stationary waiting time distribution 
[Eq. ((TB])]. The stationary weights [Eq. ((TTJ] read p^" = 0.875 
and p"^ = 0.125. The noisy curve corresponds to the time av- 
erage Eq. pO[). The parameters are the same than in Fig. 6. 



where WRit) is defined by Eq. (|63[) after the replacement 
IR -^ IR [Eq- (|S5)) ]. While Eq. ([55]) is an approximation 
valid in the limit of slow environmental fluctuations, here 
Eq. ([75]) is valid independently of the values of any of the 
parameters that define the system evolution, Eq. (|67p . 

In Fig. 8(a), we plot Wst{t,t') (as a function of t — 
t') for different values of the parameters p%{t'). As in 
the previous case, Wst{t,t') has a strong dependence on 
the values of the configurational populations {p%{t')}, 
implying strong departures from a renewal process. In 
fact, notice that depending on p%{t'), Wst{t,t') may or 
not to develops oscillatory behaviors. 



4- Stationary waiting time distributions 



The first stationary waiting time distribution Eq. (jlOp 



from Eq. 

Eflff IR (^ 



can 



-RR' 



be 

\Pr] 



written as w^ (t) = 
\+) . After using the defi- 



w 



(1) 



(r) =^Wi?(T)p' 



R I 



(76) 



where the weights p^ are determine from the relation 
>lph|poo) = I-) (-1 J2rPr\R): delivering 

„oo _ 1r{+\Pr\+)+J:r'1RR'{+\p'S'\+) ,--^ 

In Fig. 8(b) we plot the analytical expression for 
Woo (t) [Eq. (|76p ] joint with the numerical distribution 
(noisy curve) obtained as the probability distribution of 
the time intervals between successive photon emissions 
along a single trajectory, Eq. pH)) . The theoretical and 
numerical results match between them. 

The second waiting time distribution w^ (t2 , ti ) fol- 
lows from Eq. (|33l) . After some calculations we get 



w, 



(2) 



(t-2,Ti) ^ ^<WR{T2)qR + ^WR^{T2)qR'R\ 
R [ R' J 



xwr{ti)p'^. 
Here, we introduced the factors 



_ IR 

qR = -^, qwR 

IR 



IR'R 
IR 



(78) 



(79) 



which for any R satisfy the normalization qR + 

The physical content of Eq. ([78| can be read as follows. 
After a first photon emission [contribution wr{ti)p'^\^ 
the second one happens without a configurational tran- 
sition with probability g/j, while with probability qRi r it 
is endowed with the configurational transition R ^ R' . 
This interpretation is consistent with the results pre- 
sented in Ref. [3a| (see also Appendix A). On the other 

(2) (2) 

hand, while in general Woo (t2,ti) 7^ Woo (ti,T2), here 

(2) 
for two-dimensional configurational spaces Woo (t2,ti) is 

symmetric on its arguments. 

In Fig. 9(a) we plot ^^^2,^) [Eq. dTS])] for differ- 
ent values of ri. We also show the numerical distribu- 
tion, Eq. (fTTj) . In Fig. 9(b) we plot the dimensionless 
parameter A(r2,ri) [Eq. p^ ] determine from Eqs. ([7S)) 

and ([75)) . for different values of ri. In Fig. 10 we show its 

(2) 
contour plot as well as that corresponding to Woo (t2, ti). 

For small values of of T2 and ti, A(r2,Ti) is almost null, 

while for higher values of both times A(r2, ri) reaches its 

maximal values. 

(2) 
We notice that the structure of Woo (t2, n) is very sim- 
ilar to that shown in Fig. [5j The same affirmation is 
vaHd for the corresponding Woo (t), i.e.. Fig. ^h) and 
[5{b). Nevertheless, here A(r2,ri) [Fig. [TO] develops a 
much richer structure or dependence in T2 and ti. The 
origin of this characteristic can be related with the behav- 
ior of the underlying photon-to-photon emission process. 
In fact, in the light assisted case Wst{t,t') [Fig. Hl^a)], 
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FIG. 9: (a) Stationary two-time waiting time distribution 
wio (t2,ti) [Eq. (|78p ] for different values of n. From top to 
bottom we take firi = 2, 0.85, 5, 0.25, and 10. The noisy 
curves correspond to the time average Eq. (Illf) . (b) Param- 
eter A(r2,ri) [Eq. H^l] determine from Eqs. ([73 and dTSIl . 
for different values of n. From top to bottom Q.t\ — 10, 0.25, 
0.85, and 2. For clarity, the curves corresponding to the last 
three values where shifted by —0.25, —0.75, and —1 respec- 
tively. The parameters are the same than in Fig. 6. 



depending on the values of the parameter p^(i'), may 
develops strong oscillatory behaviors, while in Fig. ^a.) 
the behaviors are almost monotonous. Independently of 
the underlying environmental dynamic, by increasing the 
external laser intensity the renewal departure measure 
A(r2,ri) develops a richer structure. 

When the radiation pattern develops a blinking phe- 
nomenon '31|. i.e., for slow (light- assisted) environment 
fluctuations, {jb.R'} ^ {ir}7 Eq. ([77| becomes 



Pr 



IrPr 



Er'Ih'PI?' 



(80) 



where 1^ follows from Eq. ([50]) under the replacement 
JR — ^ IR [Eq. ([Mt ]. and the probabilities P^ are the sta- 
tionary solution of a classical master equation obtained 
from Eq. (j4ip under the replacement (^rr' — > ^r'r = 
Qr'rIr (see Eqs. (81) and (82) in Ref. ^). In the limit 
of fast environment fluctuations, {^rr'} ^ {7fl}: the 
photon emission process becomes a renewal one. 




10 15 20 25 30 35 



FIG. 10: Contour plot of wL (t"2, ti)/^'^ [upper panel] and 
A(r2,ri) [lower panel] corresponding to the plots of Fig. 9. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we formulated a quantum-jump ap- 
proach for describing the radiation patterns of single flu- 
orescent systems coupled to complex fluctuating environ- 
ments. Our results rely on a density matrix formulation 
of the problem. The master Eq. ([3]) take into account 
both the system dynamic as well as a width class of en- 
vironment fluctuations. 

The quantum-jump approach relies on a quantum mea- 
surement theory. Here, after introducing general mea- 
surement transformations acting on the system and the 
conflgurational bath space [Eq. (|13p ]. the density matrix 
evolution was written as an average over measurement 
trajectories, Eq. ((20|) . The weight of each trajectory is 
measured by its associated n-joint probability, Eq. (j22l) . 
The hierarchy of these objects completely characterizes 
the statistical properties of the measurement processes. 
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Its functional form in an asymptotic time regime pro- 
vides information about observables defined from a time 
average along a single measurement trajectory. Eqs. (|32|) 
define the stationary probability distribution for the time 
interval between consecutive measurement events, while 
Eq. p3|) define the joint probability for two consecutive 
intervals. These two objects allow measuring the depar- 
ture of the measurement process from a renewal one. 

The decomposition of the density matrix evolution into 
a set of measurement trajectories leads to a stochastic 
representation of the system dynamics, Eq. ([M)) . Each 
stochastic realization can be related with a particular 
measurement trajectory. Their structure allowed us to 
define how the measurement process occurs event-to- 
event. The waiting time distribution Eq. (j39p defines 
the probability density for consecutive recording events. 
In contrast with a renewal process, it depends on the 
stochastic state of the system, property that breaks the 
renewal character of the measurement process. This de- 
pendence encodes the influence of the bath fluctuations. 

The case when there exist only one measurement pro- 
cess, providing information about the photon emission 
events, was analyzed in detail. Independently of the 
underlying bath dynamics (and the measurement pro- 
cesses, see Appendix A) the photon-to-photon emission 
process is defined by a stochastic waiting time distribu- 
tion that parametrically depends on the configurational 
bath populations. The analysis based on Eq. PO]) allows 
to describing situations like spectral diffusion process, 
lifetime fiuctuations and molecules diffusing in a solu- 
tion. The stochastic waiting time distribution, Eq. (|5^ . 
and the first and second stationary waiting time distri- 
butions, Eqs. (|53p and ([55]) respectively, provides a deep 
characterization of the photon emission process. These 
general expressions assume a simple form when the en- 
vironment fluctuations are much slower than the optical 
system transitions. In fact, in such a case those objects 
can be written as linear combinations of the waiting time 
distribution associated to a Markovian fluorescent sys- 
tem characterized by the parameters corresponding to 
each configurational state, Eqs. ([59)) . (|64| . and (|65|) . The 
case of light assisted process, Eq. (|57|. admits a similar 
description, Eqs. dTS]), ^, and (1751) . 

The developed results provide an alternative theoreti- 
cal tool for analyzing single fluorescent systems coupled 
to classically fluctuating environments. In fact, the ex- 
plicit analytical characterization of statistical observables 
like the stationary waiting time distributions, Eqs. (jlOp 
and pTj). and the renewal departure function, Eq. (IT^ . 
may provide a power tool for deducing the underly- 
ing structure of complex nanoscopic reservoirs analyzed 
through fluorescence spectroscopy. 
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Appendix A: Measuring photon emissions and 
configurational transitions 

In Section IV we characterized the quantum jump ap- 
proach (for both self-environment fluctuations and light 
assisted processes) when the measurement action only 
gives information about the photon emission events. 
While that is the standard situation in SMS, the formal- 
ism developed in Section III allow us to analyze the case 
in which there exist extra measurement channels (ap- 
paratus) that give information about the configurational 
states of the reservoir. Besides its theoretical interest and 
potential applications, the following analysis also allows 
to understand some previous results 36]. 

Here, we assume that at any time one know which is 
the configurational state of the bath. Under this condi- 
tion, the stochastic dynamic of \pf) and \Pt^) assume the 
structure 



|pf)=p|*(i)|i?f), 



|Pf ) = \Rf), 



(Al) 



where Trs[pg(i)] — 1. and Rf randomly change over 
the set of possible values R = 1, 2, • • • i?max- Therefore, 
here the vectorial nature of \pf) can be avoided. In 
fact, all relevant information is encoded in p^git) and 
(i?|Pf ) = (5flfl.=t [see Eq. §7}]. While the underlying 
master equations are different, the results of Ref. [33| 
also rely on the previous assumption. 



1. Self-fluctuating environments 

First we analyze the case of self-fluctuating environ- 
ments, Eq. (j40p . The parameter p, includes one term 
corresponding to the photon detector, p = ph, and 
/i — 1 • • • -Rmax terms that detect (measure) when a tran- 
sition to a given conformational state R happens. 



a. Measurement operators 
The measurement operators [Eq. (jl3[) ] read 



■A^phjp) 
Mr\p) 



Hr'Ib.' T^rsiaUpii'V 
\R)J2r'^RR'PR' 



Er" 



Bf. VRR" 



Trs[pi: 



(A2a) 
(A2b) 



where \p) — ^^ \R)pii. The (unnormalized) conditional 
evolution [Eq. (|44p ] is diagonal in the i?-space and reads 

^^ = ^[HR,p\{t)] - ^R{D,p\{t)}+ - 4>Rp'k{t), 

(A3) 
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where the rate 0^ is defined by 



1.04 



(A4) 



These definitions provide a splitting of Eq. (HUl) that al- 
lows to formulate the quantum-jump approach, Eq. (1141) . 
A^ph corresponds to the transformation associated to a 
photon detection event. On the other hand, Mr take in 
account all transitions R' ^ R that leave the bath in the 
configurational state R. 



b. Stochastic dynamics 



The measurement operators Eq. 
formations [\pf) = Ps{t)\Rf)] 



2]) imply the trans- 



\pT) 
\pf) 



-A^ph|pf)-|-)(-||i?f), (A5a) 

MR\pf)=pf{t)\R). (A5b) 



The first transformation collapse the system to its ground 
state and does not affect the configurational state. The 
measurement operator Mr leaves invariant the system 
state Pg (Oi while produces the configurational transition 
\Rf) — )■ \R). On the other hand, notice that the dynamics 
between recording events, i.e., Eq. (jASp . does not affect 
the configurational bath state. 

From Eqs. HH) and (|M1) . it follows 



(A6a) 
(A6b) 



Fphilpf)] = 7flf (+|pf(t)l+): 

FrUpT)] 



'JRRT 



Consistently with the classical evolution Eq. (PT|) . the 
probability by unit of time for observing the configura- 



tional transition \Rf) 



\K) 



'■'RRf 



is independent 



of the state of the system p^' (i) . 

When a recording event happens, each transformation 
[Eqs. (lASp ] must to be selected in agreement with the 



transition probabilities tp(i), Eq. (l36t . They read 



tph(t) 



ti^(t) 



lRf{+\p's{i)\+) 

i'RRf 

7flf (+|PsWI+) + ^flr' 



(A7a) 



(A7b) 



where (t)Rst follows from Eq. (jA4p . Notice that when a 
configurational transition happens, \Rf) — > \R)^ the dif- 
ferent possible final states \R) are chosen with proba- 
bilities t^^^=t = tR{t)IY,p,,tR,{t) = (l)RRf/4>Rf. This 
result can straightforwardly be read from the classical 
master equation ({41]) . 

In Fig. 11 we show the realizations associated to the 
measurement transformations Eq. (jASP and the evolution 
Eq. (j40p . They were build up by using the finite time step 
algorithm (Appendix D). 

Fig. 11(a) shows a realization of the upper population 
of the system {+\ps{t) \+) [Eq. ^7^]. In contrast with 
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FIG. 11: Stochastic realizations of a fluorophore system de- 
fined by the evolution Eq. (|40|) and the measurement opera- 
tors Eq. (|A5|l . (a)-(c) Realization of the of the upper popula- 
tion of the system {-|-| Ps (i) |-l-) . (b) Realization of the of the 
configurational population of the bath {R\Pf^), for R = A. 
(d) Intensity realization. In (a) and (b), the parameters are 
Qr = n, SujR = 0, 7A/n = 1.5, 7s/n = 3, <^As/n = 0.03, 
(J>ba/^ = 0.05, and uil = i^o- In (c) and (d), the parameters 
are the same than in Fig. 1. 



Fig. 1, here each event may corresponds to a photon de- 
tection event, (-|-| pg (i) |-l-) — > [Eq. (jASaP J. or to a con- 
figurational transition (vertical dotted lines), Eq. (jA5bl) . 
In these last events the upper population remains unaf- 
fected. In Fig. 11(b), we show the realization of {R\Pi^), 
for i? = A. In contrast with Fig. 1, here at all times we 
know with total certainty \{R\P^^) = 1 or 0] the configu- 
rational state of the bath. 

In Fig. 11(c) and (d) we show the realization of 
(+1 Ps (^) 1+) ^^d the scattered intensity I{t). The param- 
eters are the same than in Fig. 1. In (d), the telegraphic 
signal correspond to J2r lRiR\Pt)i where the intensities 
{Ir} are defined by Eq. ([50)) . This function assume the 
value Ir when the bath is in the configurational state R. 
The plot shows the direct correlation between the value 
of the intensity and the configurational bath state. 



18 



c. Recording process 



2. Light assisted processes 



Given that a recording event happens in the /i-detector 
at time t' , the waiting time distribution for the next event 
at time t is given by Wst{t,t', fx), Eq. (15^]) . The event at 
time t is selected with probabihties t^(i), Eq. (jA7p . Here, 
Wst{t,t' , fi) can be expressed in a shorter way through 
its associated survival probability, i.e., Wst{t,t' , fi) = 
-{d/dt)Po[t,t';Mf,\pp)], Eq. §EI- From Eqs. (HH and 

^51). we get 



Po[t,t';MpiM')]=e ^"''^* *V«=.[t-t';|-)(-|], 



and for i? = 1, 2, • • • , i?n 



Po[t,t';MR\pn^e 



-4>R{t-t') 



(A8) 



W^fl[t-t';p|*(i')]. (A9) 



Here, Wii[t; p] = Tr5[p'^(t)], where p^{t) is the solu- 
tion of the equation {d/dt)p\{t) = —{i/h)[Hji,pj^{t)] — 
Jr{D, p\{t)}+, solved with the initial condition /5'^(0) = 
p. Thus, W^fl[t;p] is the photon survival probability of a 
Markovian system that begins in the state p, and whose 
characteristic parameters are 7^, ujr, and flu. In fact, the 
waiting time distribution Eq. (j60p can also be written as 
WR{t) = -{d/dt)WR[t;\-) {-\]. 

The interpretation of the survival probabilities 
Eqs. (|A8[) and (jA9l) is very simple. The exponential 
factors take into account the probability of not having 
any configurational transition in the time interval (t',t). 
On the other hand, the factors defined by WrII; p] mea- 
sure the probability of not having any photon emission 
in {t',t). In Eq. ^M^, WR[t - t';pf (i')] is the photon 
survival probability of a Markovian system (with param- 
eters corresponding to the configurational state R) that 
begins in the (arbitrary) state p^g{t'). Consistently, in 
Eq. (|A8P the factor 14^=* [i — i'; |— ) (— |] corresponds to 
the photon survival probability after happening a photon 
detection event at time t' , i.e., p^s{t') = |— ) (— | . Hence, 
here the associated stochastic waiting time distributions 
{wstit, t' , /i)} change when a configurational transition or 
when a photon recording event happen. Added to its de- 
pendence on the configurational state R, in contrast with 
the result of Section IV, Wst{t, t' , p) also may depends on 
the system state p%{t'), i.e., its functional form depends 
parametrically on the matrix elements of p%{t'). 

We have checked that the stochastic dynamic of |p^*) 
defined by Eqs. ([Ml) and ([HI) , like in Fig. 2, also re- 
cover the density matrix evolution defined by Eq. (|40p. 
As the dynamic of the configurational states is classical, 
the statistical properties of the photon-emission process 
remain the same. This fact is clearly seen by comparing 
Fig. 1(d) and Fig. 11(d). In both cases the intensity is 
characterized by the same telegraphic behavior. 



Here we analyze the quantum-jump approach associ- 
ated to Eq. (j67p when both the photon emissions and the 
configurational transitions are measured, Eq. (jAip . 



a. Measurement operators 
From Eq. (pT)) the measurement transformations read 
Y.r1r\R)^PR^'^ 



M^h\p) = 
Mr\p) - 



Efl'7fl'Trs[CTtcrpfl,]' 

\R)Y.R'lRR'^PR'^'' 
J2r" 7fli?/'Trs[fTtcrpfl/,] ■ 



(AlOa) 
(AlOb) 



R E (l,i?max), while the conditional evolution here is 
also defined by Eq. ((70)) . 



6. Stochastic dynamics 



The transformations Eq. (|A10p imply the transforma- 

tions [\pf)^p%\t)\Rf)] 



\pf) ^ >lph|pf) 
\pf) ^ MrIpT) 



\RT) 



\R). 



(Alia) 
(Allb) 



While A^ph collapse the system to its ground state and 
leaves invariant the configurational state, the superop- 
erators Mr produces both the system collapse and the 
configurational transition Rf — > R. Therefore, here any 
recording event (due to Mph or to Mr) implies a photon 
detection event. 

The transformations defined by Eq. (jAlip must to be 
selected in agreement with the transition probabilities 
V(i), Eq. dSSl). From Eqs. dH]) and (KTO\ . we obtain 



rA\pf)] = 7flf (+|P5'(i)H 

FrUpT)] = 7««f (+|pfW 
Then, the transition probabilities read 



tdct(i) 



7flr 

7i?=' 



iRit) 



iRRf 
iRf 



(Al2a) 
(Al2b) 

(A13) 



Notice that these objects are independent of the state 
p^s{t). Furthermore, they are stretched related with the 
definitions Eq. (|79l) . 

Fig. 12 shows the realizations associated to the mea- 
surement transformations Eq. (JAIII) and the evolution 
Eq. (j67p . The realizations were determined by using the 
finite time step algorithm (Appendix D). The parameters 
are the same than in Fig. 6. 

Fig. 12(a) shows a realization of the upper popula- 
tion of the system (-|-| pll(t) |-|-) [Eq. §7])]. The vertical 
dotted lines correspond to the time where the configu- 
rational transitions happen. In Fig. 12(b), we show the 
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FIG. 12: Stochastic realizations of a fluorophore system de- 
fined by the evolution Eq. (|67|) and the measurement trans- 
formations Eq. (|A11|1 . (a) Realization of the of the upper 
population of the system (+| Ps(t) |-l-) . (b)-(c) Realization of 
the of the configurational population of the bath {R\Pt^), for 
R — A. (d) Intensity realization. The parameters are the 
same than in Fig. 6. 



realization of the configurational population (A|Pf*). In 
contrast with Fig. 11(a) and (b), we notice that here 
the configurational transitions are always endowed with 
a photon emission. Furthermore, at all times the config- 
urational state of the bath is known with total certainty 
[{R\Pf) = 1 or 0] [compare with Fig. 6(b) and (c)]. 

Fig. 12(c) shows {A\P^^) over a larger time scale while 
in Fig. 12(d) we show the scattered intensity I{t). As ex- 
pected, the intensity realization is similar to that shown 
in Fig. 6(d). 



c. Recording process 

The stochastic waiting time distributions Eq. (pQ)) . 
from Eqs. 1^ and (jAlip . here read 

Wstit,t',ph) = wup{t-t'), (Al4a) 

Wst{t,t',R) = WR{t-t'), (A14b) 



where wait) follows from Eq. ((63|) after the replacement 
JR — ^ 7^, [Eq. (|69| ]. i.e., they are the waiting time dis- 
tribution of a Markovian fluorescent system with decay 
rate 7^, detuning Sfi, and Rabi frequency fir/. 

The expressions written in Eq. (|A14p only differ in 
their sub-index {Rp or R). After a recording event, 
the indexes must be chosen with probabilities (jA13|) . 
Therefore, Wstit,t',fj,) during successive photon record- 
ing events is randomly selected over the set of functions 
{w_R(t}j. This result recovers the analysis developed in 
Ref. [3a]. For the example shown in Fig. (12), the two 
functions WAit — t') and WB{t — t') can be read from 
Fig. 8(a) by taking ^^(i') = 1 and p^(^') = respec- 
tively. On the other hand, from Eqs. (|A13|) and (|A14|) . 
one can deduce that here the stationary photon waiting 
time distributions are also defined by Eqs. ([75]) and ([751) . 



Appendix B: Stationary n-joint probabilities 



The probabilities Eq. (|22|) define the ensemble statistic 
of the measurement process. They depend on the initial 
condition [po) • The statistical information that can be 
obtained from a time average along a single realization 
can be obtained from the stationary n-joint probabilities 
P^[t, {Ti}",{/ii}"]. They define the events statistics af- 
ter happening an infinite number of measurements events 
and that an infinite time elapsed since the initial condi- 
tion. 



PnrAr^rlA^^^}l] 



/■tjv 

lim lini / dijv-i ' 



p„MtAuri\{'^^}i] 



dti 


n+Ni 



The new time variables are defined as t = t — In, Ti = 
ti+AT — 1]\[. The measurement apparatus indexes are fii = 
I'i+N- By working in a Laplace domain, from Eq. (|27l) it 
is possible to obtain 



•D(t2-Ti) 

xJ^,,e^^'M\po.)]Fo.. (Bl) 



PrTir, {nViAt^rVi] = TTs[{l\e^'--^->J,^- ■ ■J.^e 



Here the measurement operator Ai is defined by Eq. pip , 
and the constant Foo reads 

roo = ^r^[|Poo)]=Tr5[(l|i|poo)], (B2) 

where F^ follows from Eq. (IT^ and we used Eq. (|31l) . 
J' = ^i^J^i- The stationary state |poo) is defined by 
Eq. dSOD- 

From Eq. (IBip . by performing the inverse calculations 
steps than in the derivation of Eq. (|27l) , it follows 



P^[rAnrlA^^^rl] - P^[r,r,,\M^APr„)] 



(B3) 



xu';,Jri,0;7W|poo)]roo- 
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The auxiliary states read |pri+i) — T(Ti+i, Ti)A^^JpT-.), 

where |p^J = r(rl,0)7W|/9oo)• 
The interpretation (and structure) of Eq. (jB3|) is simi- 
lar to that of Eq. p2l) . Nevertheless, here the factor roo 
takes into account the probability by unit of time of hav- 
ing an arbitrary detection event in the long time regime. 
The associated measurement operator is Ai, Eq. (|3T]) . 
Furthermore, in contrast to Eq. ()22|) . the first contribu- 
tion (waiting time distribution) in Eq. (jB3p is defined 
with the state M\poo), i.e., the state after and arbi- 
trary detection happening in the stationary regime. From 
Eqs. dm and (|BT]), the expressions Eqs. §^ and §^ 
follows straightforwardly after replacing ti — >■ ti, and 
T2 — >■ Ti -I- T2. In fact, the variables {ri}" of the station- 
ary waiting time distributions Woo [{7i}",{/ii}"] denotes 
the time interval between consecutive recording events. 



Appendix C: Averaging over realizations 

Here, we demonstrate that the deterministic evolution 
Eq. p^ is recovered after averaging Eq. ((34|) over real- 



izations of the Poisson processes N^ 



First, by using that {dN^ 



l^\k _ 



dNf and the property 



dN^dN]^ = S^^'dNj^, it is possible to get the relation 



E{{Nt})dNt = Ei{Nn)Trs[il\JM')] dt, (CI) 

where EdNj^}) is an arbitrary function of the Poisson 
processes {Nj^}. This equality can be immediately de- 
duced by introducing a series expansion of S. Now, we 
split the average of Eq. (pi]) as 



dt 



\pt) 



dt 



\Pt) 



d 
di 



\Pt, 



(C2) 



M 



where the first contribution is associated to the condi- 
tional deterministic dynamics and the second one with 
the disruptive measurement changes. Then, trivially it 
follows 



dt 



\pt) 



^V\p,)-Trs[{l\D\pf)]\pf). 



(C3) 



T) 



On the other hand, by using the definition Eq. ([T^ and 
the relation Eq. (JC1|) . we get 



1. Infinitesimal time step algorithm 

In the first algorithm, the stochastic state Ipj^^^j is 
obtained from |pf ), where At is the time discretization 
step. By defining the quantity 

r(t)^^r^[|pf)]=5]Trs[(l|i'^|pf)], (Dl) 



the probability AP of having a measurement event is 
defined by AP = Ai F{t). Then, a random number r in 
(0,1) is generated and compared with AP. If r > AP, no 
recording event happens, so the vectorial state evolves 
deterministically as [Eq. (IT7)) ] 



\pt^t)^nt+dt,t)\pf) 



(l + PAQIpf) 
l + Trs[(l|PAt|pf)]' 



(D2) 



If r < AP, there is measurement event. Then, the system 
state dXt + M is defined by [Eq. (IHl) ] 



\pUt)-M,\pf)- ^'^'^"^ 



(D3) 



Trsmj,\pf)] 

Here, the index p is chosen with probability t^(t), 
Eq. dSll). Due to the relation F ,,[\pt)] = t^(t)r(t), the 
generated realizations satisfy Eq. (|34)) . 



2. Finite time step algorithm 

An alternative and more efficient algorithm can be de- 
fined by using the survival probability Eq. p^ [see also 
Eq. p8|) ] . Given that the state of the system after a mea- 
surement at time ti is given by Ai^^\pf'.), the time t^+i 
of the next event is obtained from the equation 

Po[Wi,t.; \pf;)] = Trs[(l|e^(''+^-'')>l^Jp?:)] = r, 

(D4) 
where as before r is a random number in the interval 
(0,1). For t € (ii+i,ti), the stochastic state evolves de- 
terministically as [Eq. pT)) ] 



dt 



\Pt, 



M 



jApt)-Y.^^s[{i\Mpf)]\pf) ■ (C4) \pf)^f{t,t,)M,^pf:) 



e*(*-*')7W^Jpf) 



After introducing the relation Tr5[(l|2J|»)] = 
-E^Tr5[(l|J^|«)] in Eqs. ^^ and dCH, the evolution 
Eq. (fT4|) follows straightforwardly. 



Trs[(l|e^(*-*^)X^, \pf)] 

(D5) 
At time t = i^+i, an index p-i+i is chosen with probability 
{tfj_(ti+i)}, Eq. (|36|) . and then the sudden transformation 



Appendix D: Algorithms associated to the stochastic 
evolution 



Two different algorithms allow to build up the realiza- 
tions associated to the stochastic evolution Eq. (|M| . 



\ptj ^ M,^^Apt:, 



J,..APtlJ 

Trs[(l|J-p,+Jp?*^J] 



(D6) 



is applied. The first event follows from Eq. (|D4[) with 

-^MilPt*) ~^ IPo*)- The realizations generated with this 
algorithm are also consistent with the evolution Eq. (|34l) . 
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